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Abstract. We analyze the convergence to equilibrium in a family of Kac-like kinetic equations 
in multiple space dimensions. These equations describe the change of the velocity distribution 
in a spatially homogeneous gas due to binary collisions between the particles. We consider a 
general linear mechanism for the exchange of the particles' momenta, with interaction coeffi- 
cients that are random matrices with a distribution that is independent of the velocities of the 
colliding particles. Applying a synthesis of probabilistic methods and Fourier analysis, we are 
able to identify sufficient conditions for the existence and uniqueness of a stationary state, we 
characterize this stationary state as a mixture of Gaussian distributions, and we prove equili- 
bration of transient solutions under minimal hypotheses on the initial conditions. In particular, 
we are able to classify the high-energy tails of the stationary distribution, which might be of 
Pareto type. We also discuss several examples to which our theory applies, among them models 
with a non-symmetric stationary state. 



1. Introduction 

This paper is concerned with the spatiahy homogeneous Boltzmann equation 

dtfJ. + ^J.^ Q+[fi], /x(0) = ^o, (1) 

for a time-dependent probabihty distribution /i(i,dw) on the particle velocities v G Mf^ in a gas. 
Above, the collisional gain operator = Q+[fi,fj] is bi-linear and given (in weak form) by 

ijiv')Q+[fii,fi2]{dv') = l f E[i;{v')+^iv:)]^^l{dv)^l2{dv^.), (2) 

for every ip S Cfc(R'*). It accounts for the change in velocity due to binary collisions between 
particles. Taking the expectation E in ^ corresponds — morally — to integration with respect 
to the scattering angle. Our main assumptions are the following: 

(1) Collisions occur at unit frequency, i.e., the probability for a particle to collide is indepen- 
dent of its velocity. This justifies the particular form ((1} of the Boltzmann equation. 

(2) The post-coUisional velocities v' and v'^ are determined from the pre-coUisional velocities 
V and by linear rules 

v' — Lv + Rv,,, v'^ = L^,v^, + R^,v, (3) 

in which L, R, and i?* are d x d random matrices of a given distribution, that is 
independent of v and v^. For simplicity, we assume that (L, R) and (L*, R^) are identically 
distributed, but we do not assume their stochastic independence. 

(3) The model conserves the second absolute moment of fi{t; dv), which means that 

E[k'P + KP] =kP + k*p. (4) 
Thus, the total kinetic energy (or temperature) is a constant. 
Equation ((!]) with kernel constitutes one possible generalization of the celebrated Kac model 
to multiple space dimensions. The original Kac model — which is ^ and ([2]) in d = 1 dimension 
with {L,R) = (cos 0, sin 0), and (p uniformly distributed on [0,2tt) — has been introduced [53] 
as a caricature of the Boltzmann equation to investigate the propagation of chaos in particle 
systems. In multiple space dimensions. Maxwell molecules [5] are the commonly accepted canonical 
"interpolation" between the simplified Kac model and the genuine Boltzmann equation. Also in 
the Maxwell case, one uses linear rules of the type ([3]) , which is the so-called w-representation of 
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collisions, and a uniform collision frequency, just as we do here. The difference is that for Maxwell 
molecules, the probability distribution of the coefficients (L, R) typically depends on the spatial 
direction oi v^, — v. (Maxwell molecules without this dependence are a special case of the model 
considered here, see Section [Q] ') 

Our approach is not an attempt to bridge the gap between the Kac model and the physical 
Boltzmann equation in an alternative way. Rather, we introduce a novel class of binary interaction 
models, which are reminiscent of the Kac equation, but with distinguished features that makes 
this class interesting on its own right. We emphasize that there is a huge freedom in the choice 
of the matrix coefficients {L,R). The theory that we develop does not require the conservation 
of momentum or energy in binary collisions, but only the much weaker stochastic condition 
Neither do we ask for any covariance under rotations of M''. Binary particle collisions without 
momentum/ energy conservation have been used before, for instance to combine particle collisions 
and the influence of a heat bath into a single interaction mechanism [TT] , as well as in applications 
of kinetic theory to simple models for wealth distribution [161 112] . The violation of rotational sym- 
metry, on the other hand, can be thought of as resulting from a non-isotropic background, which 
induces an a priori preference for certain post-collisional directions on the interacting particles. As 
a consequence of its high degree of freedom, there is a huge variety of potential stationary states for 
the velocity distributions, the details of which depend sensitively on the model parameter. These 
stationary distributions might possess high energy tails of Pareto type, and are not necessarily 
rotationally symmetric. 

For the development of the theory, we continue in the spirit of various recent works, like [H El 
[T9l [20] , which have been concerned with the applicability of probabilistic methods to study the 
equilibration in kinetic systems. Specifically, our approach parallels in large parts the one taken 
recently in ^|3j for the one-dimensional reduction of ([T]). In one spatial dimension, when L and 
R are simply non-negative random numbers £ and r, it has been proven that the conditions 

E[f + 7-2] = 1, and (5) 
E[F + rP]<l for some p > 2, (6) 

are sufficient for existence and uniqueness of a stationary distribution /ioo of finite second moment. 
Moreover, if ([5]) and ([6]) are satisfied, then finiteness of the second moment of the initial condition 
/iQ is sufficient for (weak) convergence of the transient solution /x(t) to /ioo as f — cx). 

Here, we are going to prove an analogous statement in the multi-dimensional situation d > 1. 
However, the substitutes for conditions ([5]) and ^ are not obvious. In fact, the most straight- 
forward generalizations of these conditions to matrices L and R would not work. Moreover, our 
proof needs a third condition, see ((T2|) below. We are now going to state and briefly discuss the 
individual conditions. 

The first condition, which parallels (O, accounts for our hypothesis that the second absolute 
moment is constant in time. 

Assumption 1. The matrix coefficients L and R satisfy 

E[L'^L + R'^R] = 1, (7) 

where 1 G M.'^^'^ denotes the unit matrix. 

It is tempting to substitute ^ by the condition 

E[\\L\\P + \\R\\P] < 1 

with, say, the operator norm ||A|| — sup|g|^]^ In general, however, this condition would be 

too strict. It is not even satisfied in the paradigmatic example of Maxwell molecules discussed in 
Section [6. II Instead, the following appears to be the most appropriate substitute. 
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Assumption 2. There exist real numbers p > 2, lo > \ and Kp e (0, 1), and a weight function 
a; : M"^ — > M with the following properties: for every rj G R'*, 

a;(A7/) = XPuj{t]) for all A > 0, (8) 

|??r < ^iv) < ^M'', and (9) 

E[uj{L^7]) + ujiR^v)] < Kp ^iv)- (10) 

We refer to ^ and ^ as p-homogeneity and boundedness of uj, respectively. 

This condition is employed at various points of the discussion; in combination with Assumption 
[3] below, it provides contractivity estimates in suitable metrics and weak compactness of transient 
solutions. Notice that w(^) = defines a p-homogeneous and bounded weight function. In 
Section 16.31 we give a multi-dimensional example in which a more complicated weight function is 
needed in order to satisfy (fTO|) . 

Our last hypothesis has no analogue in one dimension. It guarantees that the multi-dimensional 
model cannot be decomposed into a family of one-dimensional sub-systems. 

Assumption 3. The map T, defined on symmetric matrices M G K.''^'^ hy 

T{M) = E[LML^ + RMR^], (11) 

possesses a fixed point = T(S*) that is positive definite and has trE* = d. Moreover, there is 
a constant k < 1 such that 

||T(M)||oo < K||Af||oo. (12) 
holds for every traceless symmetric d x d-matrix M . 

Here ||M||oo denotes the spectral radius of the matrix M; see Section [2721 for its definition. An 
immediate implication of (jl2p — in combination with ([T]) — is that the fixed points of T are scalar 
multiples of S*. This condition is needed to ensure that all "directional information" of the initial 
condition is lost in the long-time limit. 

Assumption [3] — although not very intuitive — appears to be essential for the development of 
the theory. It accounts for the fact that we are not working in the one-dimensional or radially 
symmetric situation, like in [3] or [6], respectively. Indeed, one could think of a straight-forward 
generalization of the one-dimensional/radially symmetric concepts from [3] or [B] to general solu- 
tions in R*^ by using 

L^il, R = rl, 

where i and r are non- negative random variables satisfying ([5])&(|ni). Clearly, while Assumptions 
[T] and [2] hold in this case. Assumption [3] is violated, since 

||E[LML^ + lloo = + r^]M\\^ = \\M\\^ 

for every matrix M. By the results of [3], one can associate to every one-dimensional linear 
subspace span(w) C R'' with 7^ w G M'^ a non-trivial stationary distribution //JJ^, which attracts 
all transient solutions fi{t) to (H]) that are initially supported on span(w) and have unit second 
moment. Consequently, there are infinitely many stationary distributions. 

We shall now summarize our main results. Below, we mean by temperature of a centered 
probability distribution /z on M.'^ the quantity 

The normalization is chosen such that the Gaussian distribution with unit covariance matrix is of 
unit temperature. 

Theorem 1. Under Assumptions 1-3, equation ^ possesses precisely one stationary distribution 
A*oc = Q+[Moo] that is centered and of unit temperature. This distribution fioa is a scale mixture 
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of Gaussians, i.e., 

for every test function ip G Cb{M.'^), where S is a symmetric and positive semi-definite random 
matrix m M''^'', and ^/S denotes its positive semi-definite square root. Moreover, the pth absolute 
moment of /loo is finite. 

If fiQ is a centered initial distribution with finite temperature, then the corresponding transient 
solution of the initial value problem ((IJ converges weakly to a non-trivial stationary distribu- 
tion fi^, and fi^ is obtained from fioo by the unique isotropic rescaling for which 0[/i*] = 0[/io]. 

// in addition P{L-^e — R^e = 0} = and F{L'^e ^ Q ^ R^e} > for every unit vector 
e, then /^oo is absolutely continuous. In this case, if fXo is a centered initial distribution then the 
corresponding transient solution fi(t) converges weakly to /i* if and only if has finite temperature. 

Apart from Theorem [1] we prove results on the Sobolev-regularity of /j,oo 's density (Theorem 
HI, on the rate of convergence of ^[t) to /ioo in Fourier metrics (Theorem [5]), and on the divergence 
of higher moments (Proposition [3]). 

FinaUy, we outhne the structure of the paper and give an overview on the apphed techniques. 

Section [5] coUects various definitions and prehminary observations. In particular, we introduce 
a family of matrix norms. 

Section [3] is concerned with the existence and certain properties of the stationary distribution 
/^oo- Assuming that /ioo is of the form p3p . we first show that the random matrix S needs to 
satisfy the fixed point equation 

Law(S') = La.w{LS'L^ + i?S"'i?^). (14) 

The right-hand side of ([H)) — where S' and S" are i.i.d. copies of S, which are independent of 
{L, R) — is a multi-dimensional generalization of a smoothing transformation. For (scalar) random 
variables, fixed points of such transformations have been intensively studied in the literature, see 
e-g-, [m [221 mi US] . Some multidimensional generalizations of the smoothing transformation have 
been treated in [1] [30] , although there are apparently no results on equations of the form ([T^ . 

In Theorem [21 we obtain the (essentially unique) solvability of (fH)) under Assumptions 1-3. 
The proof uses a generalization of Fourier metrics. Originally, these define a distance between 
probability measures on M.'^ from a suitably weighted difference of their respective Fourier trans- 
forms; see pi)] for a survey. The basic ideas have been developed in [5] and [H], and since then 
applied — with appropriate modifications — in various context related to homogeneous kinetic 
equation of Kac type, see e.g. 01 \W\. Our generalization to probability measures on matrices 
uses the same basic concept, but is technically more involved. 

Subsequently, results on symmetries, moments and further properties of /ioo are derived. The 
most interesting auxiliary result is probably that on the existence and regularity of a density. 
Since /too is not single Gaussian in general, but a mixture of the form (|27|) . it might possess 
concentrations. In one spatial dimension, one can generally decompose any such scale mixture 
into a singular measure concentrated at the origin, and a regular measure with a density that is 
C°° on R \ {0}, but in general only on R. We refer to [27] for a discussion on the Sobolev 
regularity of stationary states, and on the propagation of smoothness in ([1]) with ([2]), and to 
[21 [131 [TB] for closely related results on Maxwellian molecules. 

The situation is more difficult in multiple dimensions, since concentrations can occur not only 
at the origin, but on any lower-dimensional linear subspace, or on the union of subspaces. A 
non-trivial example is given in Section 16.31 We do not discuss the variants of possible concentra- 
tions in detail, but instead provide a sufficient condition for absolute continuity of the stationary 
distribution in Theorem [31 Its proof is a multi-dimensional extension of the particularly elegant 
approach made in j26[ . 

In Sectional we prove that arbitrary solutions fj,{t) to the initial value problem ([T]) converge to 
/ioo if they are initially of unit temperature. We use two approaches that lead to slightly different 
results. 



-dw 



(13) 
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The first result, see Theorem |4j is obtained by probabilistic methods. The key element is a 
stochastic representation of fj,(t) as the law of a weighted sum of independent random variables. 
The long-time asymptotics of ii(t) can then be obtained by studying the convergence of these 
weighted sums. The latter is conveniently done in the framework of the central limit theorem. 
The idea goes back essentially to McKean [55J [55], and has been brought to full power in the 
recent works |14j [15] . We are able to adapt large parts of the general strategy — that has been 
successfully employed in a related one-dimensional situation [5] — to our current needs. The 
result is qualitative as it gives no direct information on the rate of convergence of fi{t) to /Xqo; the 
advantage of this proof is that it works under the minimal hypothesis on the initial condition fj,Q, 
namely finiteness of the second moment. 

The second result on convergence in Theorem [5] makes the previous statement quantitative as it 
provides an exponential rate for the convergence of nit) to /Xoo in Fourier distance. The additional 
requirement is finiteness of the p'th absolute moment of /xq for some p' > 2, and the exponential 
rate depends sensitively on p' . The proof, which closely follows a by now classical strategy, is only 
sketched. 

The proofs are carried out for an initial datum of unit temperature, 0[/io] = 1, in which case 
the associated transient solution converges weakly towards the stationary distribution fi^o with 
covariance matrix Cov/ioo — S,, the fixed point of T. A simply scaling argument shows that 
if -d :— 0[/io] is not one but finite, then the transient solution still converges to a stationary 
distribution /i^, which is simply a rescaling of fioo, with covariance matrix Cov/xJ^ = 

In Section [5] we analyze the limiting case 0[/io] = +oo. The transient solution ^(t) exists 
for all times, but it does not converge weakly to a limit measure as t — > oo. Loosely speaking, 
fi(t) concentrates at infinity instead. This statement has been made precise in the case of the 
spatially homogeneous Boltzmann equation for pseudo Maxwellian molecules in |2 and for the 
Kac equation in jj]: the restriction of fj,{t) to any bounded set in phase space vanishes in the long 
time limit. 

In Theorem [6] we prove — based on the ideas from [Zl [8] — that the situation is the same 
for the class of equations considered here. The main idea is to consider a restriction fiQ of the 
initial condition /iq to a large ball, such that /Iq has finite but huge temperature t9 — &[p.i] ^ 1, 
and compare the original transient solution to the one associated to the initial value fLQ. 
Combining the probabilistic representation with the quantitative convergence results from Section 
[31 we thus obtain an upper bound on the density of fi(t) — for sufficiently large t — by a multiple 
of the rescaling /ij^ of /ioo- For — ^ cxj, this upper bound vanishes, and consequently, the density 
of must converge locally uniformly to zero as t — > oo. 

Finally, Section [B] provides several examples to which the developed theory applies. 



2.1. Conventions. For simplicity, we shall assume throughout the paper that the random co- 
efficients L and R are symmetric in the sense that {L,R) = {R,L). Above and in the rest of 
the paper, Z = Z' means that the random elements Z and Z' have the same distribution, i.e. 
Law(Z) = Law(Z'). 

2.2. Spaces of probability measures. By ViX), we denote the space of probability measures 
on X. Introduce the class Ci{M.'^) of centered probability measures on R"*, 



of measures with unit temperature and finite qth absolute moment. The natural topology on 
C<j(R'^) is the one induced by weak convergence of probability measures, plus convergence of the 
qth absolute moment. 



2. Notations and preliminary results 




and for g > 2 the subclasses 
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As usual, we denote by the R'^^'^ the set of aU real d x d-matrices. Introduce the subsets of 
symmetric and positive semi-definite matrices by 

5 = {S € R'^'^'^IS'^ = E}, 5+ = {S e S\v^J:v > for aU v e R'^}, 

and define accordingly for p > 2 

Cp{S+) = i^i^ eV{S+) J (tr E) i.(dE) = d, y (trE)5 i.(dS) < ooj , 

the probability measures on 5+ with finite pth absolute moment and normalized first moment. 

2.3. Matrix norms. Several different norms for matrices will be used in the sequel. First recall 
that a real symmetric matrix E G 5 always possesses a spectral decomposition of the form 

d 

i=i 

where {ej}j^^ is an orthonormal basis of of eigenvectors of S, and the Xj are the corresponding 
(real) eigenvalues. For q > 1, introduce the qth matrix norm by 

d 



as well as 



E||oo '■— max |Ai| — sup |e Ee| 
1=1,. ...d 



All norms |1 ■ IL are equivalent. In particular, the 2-norm is induced by the scalar product 



9 

d 



tr(S^E) = J2 ^^i^^i- 

The following submultiplicativity estimates will be used, 

|tr(S^E)|<||S|li|lE|U, (15) 
\tv{A^^A)\<\\n\2\\A\\l (16) 
\v^^w\<m\Mwl (17) 
which hold for arbitrary E, S G 5, A G M''^'^ and w, w G W^. 

2.4. Weight functions. In several proofs, we shall work under the simplifying hypothesis that 
Assumption [2] holds with p G (2, 3). We shall now show that there is no loss of generality in doing 
so. 

Lemma 1. For every p' G (2,p), Assumption\^ is satisfied with the modified weight function 

Oj'iO^C'ojiO'-'ie^^O', where e:=^G(0,l), 

p-2 



an 



d constants lo' — {C^j^ / C,-Yu)^ ' ^ 1; 

Kip' — Ki p ' G (0,1), where C+iC- ^ ^+ '""e ^he biggest and 
smallest eigenvalue ofY,<,, respectively. 

Proof. One has a;'(AO = X^^^'^^ X^^uj' [C) = Ap'w'(0 for aU ^ G M'^ and A > 0, showing ©. And 
from C-ieP < e^*^ < C+\e it follows that (C-fC-m''' < ^'(0 < (C+/C-)^<:^'-'Cier', which 
is (O. Finally, by Holder's inequality and the symmetry {L,R) ~ {R,L), 

E[oj'{L'^0+^'{R^O] = 2E[w'(L^0] < C'(2E[c^(i^e)])'"'(2EK^iE,L^e])' 
This confirms ^Wh. □ 
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From the weight function w on R, define the following weight function on symmetric matrices 

d d 

= ■mi[j2^{^j) I S = ^r,e,ef with R'', n, . . . ,r, € {±1}}. (18) 

The next lemma lists useful properties of W; recall the definition of p > 2 from Assumption [5] 

Lemma 2. The weight function W 

(i) is homogeneous of degree p/2: 

W(AS) = AP/2>V(S) /or S e 5, A > 0; (19) 

(ii) is compatible with the usual matrix norms on the set S: 

^-(p/2-l)p||P/2 <_p||P/2 fg^all^^^. (20) 

(iii) satisfies the following analogue of (jlOp ; 

E[W{L^EL) + W{R^ER)] < KpW(S) for all EeS. (21) 
Proof. Relation (|19p is an immediate consequence of the homogeneity of w and the fact that 

n TL 

To prove the second inequality in (j20p . let S G 5 be given with spectral decomposition 

d 

2 = (22) 

Choosing = -^/jAjfej and tj — sgn Aj , it follows by definition of W as infimum and in combina- 
tion with © that 

d d / d \ p/2 

Wis) < ^u;(o) <zjj2 lA.r^' < ^ E i^^i 

To obtain the first inequality in ((20)) . choose S € 5 and e > 0, and let ^i, . . . , € R'' be such that 

d d 

2 = E^^-0ff and ^c.(e,)<W(S) + e. (23) 
In combination with the spectral decomposition (P^. one has 

iisiii = X:i^^-i-Ei-^--^-i = E E-.(-i^o)' < E(erc.)^ = Eic.p. 

fc=i fc=i fc=i j=i i,fe=i i=i 

Then, using Jensen's inequality — recall that p/2 > 1 — and the boundedness Q, 

Taking the limit e | shows (PO)) . Furthermore, again under hypothesis (1^^ . one has for every 
matrix pair (L, R) that 

L^EL = ^r,(L^O)(i^Gr, R^^R = E-.(^'^0)(^^e,r. 
And consequently, 

d d 

E[W{L^EL)+W{R'^ER)] < ^E[uj{L'^S,j)+oj{R'^£,j)] < Kp^uj{^j) = KpW{E) + Kpe, 
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proving ((2T|) in the limit e | 0. □ 

3. The stationary problem 

In this section, we prove — under Assumptions 1 to 3 — the existence of a non-degenerate 
distribution fi^o € Cp{R'^) satisfying the stationary Bohzmann-equation jioo ~ Q+[fJ'oo]- Also, 
we analyze the regularity of iioo and investigate the finiteness of its moments above the second. 
The existence theory is based on the ansatz ([25)1 below; we defer the proof of uniqueness for the 
stationary solution to Section H) 

3.1. Existence of a stationary distribution. Throughout this section, we work under the 
hypothesis that Assumption [5] holds with some p e (2,3); by Lemma [T] above, there is no loss of 
generality in doing so. Next, recall that the Fourier- Stieltjes transform /i of a probability measure 
on M'' is defined by 

m - / e^«"V(d«), 

for every ^ g M''. Taking ip{v) — exp(i^ • u) as a test function in ([2]), we obtain the following 
Fourier representation of the gain operator: 

g+[Ai,A2] = E[Ai(£^0A2(i?^0]- (24) 

We prove the existence of a solution ^/ — /too to the stationary equation 

vI/(0=E[vIt(L^O*(^''0] (25) 
which is a scale mixture of Gaussians, 

*(0=E[exp(-i^^5^)]. (26) 
Here 5 is a random matrix in 5+, whose distribution is to be determined. 

Remark 1. The representation (|26p is equivalent to (|13p in Theorem QJ Assuming that S is 
positive definite a.s., then (I13p takes the probably more familiar form 

d/Zo 



du 



-(v) =E 



1 / 1 T -1 

, exp V S 1 

det V2^ V 2 



(27) 



of a (non-isotropic) scale mixture of Gaussians. 



Lemma 3. Suppose that i>oo is a probability measure on 5+ with the following property: if S , S" 
and S" are random matrices with identical distribution Voo, which are mutually independent and 
also independent of{L,R), then 

S ^ LS'L^ + RS"R^. (28) 



Then (j26p defines a solution to 

Proof. Let B be the cr-algebra generated by {L,R). Define ^E" as in (I^Bl) . Then, using 

E[*(i^O*(^^0] = E[E[exp ( - ^{L'^0'^S'{L^O)\S]E[exp ( - ^{R^O^ S"{R^O) \B] 

= E[cxp ( - \(^{LS'L^ + RS"R^)0] 
= E[exp(-ie^5C)] =vi/(^) 

holds for every ^ G M'*. □ 

The rest of this section is devoted to proving 

Theorem 2. Under Assumptions 1-3, there exists precisely one measure v^o G Cp{S) which 
satisfies the hypotheses of Lemma\^ Consequently, the probability measure jirx, onW^ associated 
to t'oo via (|26p is a solution to the stationary equation = ^i in the class Cp{W^) . In particular, 

/ioo has unit temperature and finite pth moment. 
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We shall obtain i/qo as the fixed point of a map T defined as follows: Given two probability 
measures v' ,v" on iS+, let S' and S" be random matrices in iS-j. with distribution v' and v" , 
respectively, that are independent of each other, and independent of {L,R). Define 

T[iy\ v"] = Liiw[LS'L^ + RS"R^]. 

Clearly, any solution Voo of = T[j/, i'] verifies the hypotheses of Lemma [3l 

Lemma 4. Let Assumptions 1 and 2 hold. IJ v' and v" belong to Cp{S-^-), then T[i>' , v"] £ Cp(iS-f). 

Proof. By linearity and cyclicity of the trace, by Assumption 1, and since {L, R) = {R, L), 
trE[LS"L^ + RS"R^] = E[tr(LS"L^ + RS"R^)] = tr(E[L'^L]E[S"] + E[i?^i?]E[S"']) 
= itr(E[5' + y]) = 

Non- negativity is preserved, since S'v > for all v e K'' clearly implies that nF {LS' L'^)w = 
{L'^ wY S' {L'^ w) > for every w G R'^, and similarly for S" . Finally, note that E[ti{S')P/'^] < +oo 
by hypothesis, that E[||L^||2] < +oo follows from Assumption [5J and that hence 

E[(tr(i5'L^))^/'] <E[(tr5')'^>[||i^||^] < +oo 

by inequality (jl6p . and by independence of S' and L. The respective estimate on ti-{RS" R^) 
follows in the same way, thus concluding the proof. □ 

In order to prove existence and uniqueness of a solution to = r[i', v], we are going to show 
that T is contractive with respect to the following distance between elements i^i, G Cp(iS+): 

T>w{'^i,'^2) ■■= Sw{>'i,'^2) + a\\E[Si - S2]\\^. (29) 

Here a > is a parameter (to be specified later), 81,82 are independent random matrices with 
respective distributions 1^1 , 1^2 , and 

Sw{iyi.iy2):= sup W(S)-i|i>i(S)-i>2(S)-itr(S^E[5i-52])|. 
a^Ees 

j),;(S) denotes the characteristic functions of I'i at S G K'^^'', which can conveniently be written as 
follows: 

i>,(S) =E[exp(-itr(S^5,)], (30) 

with 8i being a random matrix of distribution i>i. In fact, it suffices to evaluate j^i^S) on matrices 
S G 5, since if S-^ = — S, it follows by the cyclicity of the trace and from 8^ — 8i a.s. that 

-tr(H5,) - tr(S^50 = tr((5f S)^) = tr(5,S) - tr(S50, 

and hence equals zero. 

Lemma 5. The distance Dyv defines a metric onCp{S+) with the following convexity property: If 
vi = sv'i + (1 — s)v'{ and V2 ~ SV2 + (1 — s)v2 are two convex combinations of measures in Cp{W^), 
with s G [0, 1], then 

'Dwiyi,v2) < s'Dw{'^'i,'^'2) + (1 - s)D>v(z^i,i^2)- (31) 
Moreover, Dyy provides a pointwise control on the Fourier transforms, 

|j>i(S)-i>2(S)| <i^max(||S||i,||S||?/')Dw(^^i,^^2), (32) 

with a constant K that only depends on a and W. Finally, convergence Dyv(z^fc,zy,) — > for a 
sequence in Cp{S+) implies weak convergence of Vk to and convergence of the expectation 
values. 
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Proof. The only difficulty here is to prove well-definiteness of Dyy: given z^i, 1^2 G Cp(5+), we need 
to show that 6w{i'i,i'2) < 00. Using estimate (llOip from the appendix in the definition ([50)) of 
i/i, we obtain (recall that we assume p G (2, 3) throughout this section) 

|i>i(S) - (l-iS^E[5i])| <E[|exp(iS^5i)-(l + iS^5i)|] 

<C,E[|S^5in <C,||S||?%[||5ir^^]. 

And thus, by the triangle inequality and (j20p . 

< CpE[\\s^rJ,^ + \\S2\\'^'] sup (w(s)-i||s||5'/2) < Cpd'^r2-i^(^tTS^r/' + {tTS^r^']. 

H#0 

The last expression is finite by the definition of Cp(iS+). 

The rest of the proof is straight-forward: Symmetry, non-negativity and the triangle inequality 
follow directly from the definition of Dyy. The convexity property pip is a consequence of the 
respective individual convexities of 6w and || • ||oo- The estimate p2p is obtained by observing 
that 

|:)i(S) - :>2(S)| < W(S)Dw (1^1,1^2) + I tr(S^E[5i - 52])| 
follows directly from the definition of Dyy, that 

I tr(S^E[5i - 52])| < ||S||i||E[5i - 5 - 2]||oo < -||S||iDw(i^i, 1^2) 

a 

holds by relation ([T5|) for matrix norms, and by recalling relation ([20)1 between and W(E!). 

Finally, if Dw(i^fe,z^*) — ^ 0, then in particular E,[Sk] E[S'h.], and the pointwise convergence 
j>fe(S) — > z>oo(2) for every S e 5 due to ([5^ implies weak convergence of the measures i^k to 
ly*. □ 



Lemma 6. Let Assumptions 1-3 hold. Let {iy'i,i>i) and (i^2'^2) ^'^'^ pairs of measures, each of 
which belongs to Cp(iS+), aKC? define vi = i^"], V2 — r[i^2, i/j]- T/ien 

5w{yi.y2) < fiplliSwii^'i, ^2) + Swi^l ^2 )) + d^/'(||IE[5; - 5^]||oo + ||E[5r - 5^']||oo) | • (33) 



where S[, 5" respectively, are independent random matrices with law v[, v'l respectively, for i 
1,2. 

Proof. Fix some E E S. For i^i — T[vl, i^"], we obtain the representations 
v^{~) = E[exp ( - Hx'EF {LS[L'^ + RS'^R^))] 



E 



E[exp ( - itv{L^ELf S[) |6]E[exp ( - itr(i?^Si?)^5f)|6] 



= E[A:(L^SL)i>,"(i?^Si?)]. 



Now define 



A' :=E[5; -5^], A" :=E[5;'-5^'], 

and observe that 
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where Si has law i^i. By the triangle inequality, using that characteristic functions are bounded 
in modulus by one, 

d:= |i>i(S)-i>2(S)-ztr(S^A)| 

= \E[{^[iL^ELy;{R^ER) - z>^(i^Si)j>^'(i?^Si?) - ztrS^(iA'L^ + RA"R^)]\ 

< E[|i>; (L^SL) - ly'^iL^EL) ^itriL^ELfA'l] + E[|l - i>;'(i?^Si?)|| tr(L^SL)^A'|] 

+ E[|i>;'(i?^Si?) - 0'i{R^ER) - ztr(i?^Si?)^A"|] + E[|l - 0'^{L^EL)\\ tr(i?^Si?)^A"|] 

<E[>V(L-.L)]sup^l^^("')-^^("')-^^^("')"^'l' 

H'#0 



sup 
e:"#o 



W(S') 

\v'{{'E.") - v'^{E") - ztr(S")^A'^| 
W(S") 

' l-v'{{E")\\ 



+ E\\\R^ml'^''W^m^^'\] sup i 
+ E\\\L'EL\\\'^-'\t,{R^ERf A"\] sup 

H'#0 



|2"||P/2-l J 

, lis'll?/-^ 



The first two terms are estimated further using E[yV(L^Si)] = E[yy(i?^Si?)] < (kj,/2)W(S), 

_____ 

which follows from (PT|) and the symmetry {L,R) — {R,L). For estimation of the last two terms, 
we apply the rule ([T5|) to the expressions containing the trace, and then Young's inequality, see 
([M)) . to the product inside the expectation. Altogether, this yields 



6< 



sup 

|P/2 



|z>J(S')-i>^(S')+ar(S')^A'| 



sup 



±E[||L^SL||r + ||i?^Si?||?/^] 



lA'l 



sup 

H"#0 



|i>;^(S")-i>^^(S^O+ztr(S")^A"| 

'ii-^2(s')r 



f ||A" 

1 / " ^" \ IM 111 



|=//||P/2-l 



oo sup 
H'#0 



"/||P/2-l 



The second terms needs more estimates. Observe that by ([20)1 and (pTjl . one has 

It remains to bound the supremum involving |1 — i'2(S')|. Since 1^2 € Cp(5+), one has trE[S'2] = d. 
Invoking (fT5)l . and using that ||S||oo < tr(S) for every S G 5+, one finds that 

|l-i>^(S)| i'^^CsS) d.s<^V[-*MS^^2)exp(-istr(S^^^))]|ds 

< ||S||iE[||5^||oo] < ||S||iE[tr^^] 
In combination with |1 — i>2(S')| < 2 for all S' G 5, one derives 

|l-'^2(-)l / :„ rol|'="||l-P/2 ,|,^/||2-p/2^ 



sup 



^ < sup min{2||S'||}-''/',d||S'||^''/'} = 22-^/2^^/2-1 < 3^, 



and the analogous estimate for the term involving j>"(S"). Putting everything together, l\'3'3\i 
follows. Note that we have implicitly used our additional assumption 2 < p < 3. □ 

Lemma 7. Under the same assumptions and with the same notations as in Lemma\^ it follows 
that 



with the constant 

A — max (^p, K + 2a-^(i^/^Kp) , 
which is smaller than one if a is large enough. 



(34) 
(35) 
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Proof. It suffices to combine equation (l33l) above with the observation that 
||E[5i - S2]\\ = \\E[L{S[ - S!,)L^ + RiS'l - S'^)R^]\\ 



<-{\\T{E[S[-S',])\ 
<'^mS'i-S'2]\\oo- 



+ \\T{E[S';-S'^])\ 



Here we have used the triangle inequahty for 1 1 • 1 1 oo , the symmetry {L,R) = {R,L), and Assumption 
[3] in combination with the fact that trE[S'J] = trE[S'2] = d. □ 

Lemma 8. Let Assumptions 1-3 hold. There exists a constant C such that the following is 
true: For arbitrary v' G Cp{S+), and if S , S' and S" independent random matrices with law, 
v:=T[v'^i'"], v' and v" , respectively, then 



sup 



< I sup ■ 



sup 



c. 



(36) 



Proof. The starting point is an apphcation of inequahty (|100p : 



sup = sup 



E 



:[{e{LS'LT + RS"R^)^Y^''\B] 



^(0 '^(O 
< sup^lE[E[ifLS'L'^0''^^\B]+E[{i^RS"R^O''^^\B]] 

+ CpE [E[{eLS'L'^m^RS"R^O''^'-' + LS' L''0''^'-\e RS" R^0\B] ] } ■ 
The first two terms inside the expectation are estimated directly using Assumption [2j 
sup ^(EfE[(^^L5'L^0P/2|S] + E[{ fRS"R:^0''^^\B]] \ 

e^n UJit) I -'J 



< supE 



< 



Kp C 



sup 



E[(e'^S"e')^^^ 



E[(e"^S"'e")P/^] 



}■ 



To estimate the remainder terms, we apply Young's inequality (|99p and inequality (jl7p . 
E 



E[(e^L5'L^C)(C^^^"i?^ef/'"'|S] <E[||5'||oo]E[||5"||L^'-^]E[|L^eni?'^er'] 
<E[tr5']E[(tr5"f/'"'] 



(37) 



p p 

Finally, Jensen's inequality gives E[{iT S"y/'^-'^] < E[tr ^ dP/2-i; recall our additional 

assumption that 2 < p < 3. A similar estimate is derived for the other remainder term; in 
summary, we obtain 

E[E[ie LS'L'^0{C^RS"R^O''^'-' + (e LS' L^O'^^-\e RS" R^0\B] 

< dP^^KpE[uj{L'^^)+uj{R^£_)] < dP^^KpUj{£_). 

This means that dM]) holds with C = d^/'^npCp. □ 

Proof of Theorem\^ Choose some z/q € Cp(5+), and define inductively the iterates Vm — T[vm-i, Vm-i] 
for all m > 1. Denote by S'o, 5*1, 82,... mutually independent random matrices with respective 
distributions 1^0, ^i, 1^2, ■■ ■ Iteration of estimate (1551) gives 

/2I 



C/o w(0 ^ 55^0 w(0 



1 - K. 



1 Kn 
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for every m > 1. This implies the m- uniform bound 

n\Sm& < ^E[||5o||g/'] + (38) 

L — Kp 

which means that the Vm form a tight family of probability measures on 5+ . Consequently, there 
exists a subsequence Vm> that converges weakly to a limit Voo G Cp(5+). 

The next step is to identify the limit Vao as a solution to the problem (j28p . i.e., we need to 
prove that Voo = T[vao, J^oo]- This is easy on the level of the Fourier transforms we know that 

i>m'+i (S) = E[i>„, (L^SL)i>™, (S e S) (39) 

for all m > 1, and we need to conclude that 

z>oo(S) = E[j)oo(i^SL);>oo(i?^Si?)] (S G 5). (40) 

Weak converges implies in particular that the Fourier transforms converge pointwise to the 
Fourier transform i>ca of the limit. Passing to the limit m' — > oo on the right-hand side of ppp . 
we obtain — using dominated convergence — the right-hand side of (I40p . for every S € 5. 
On the other hand, iteration of estimate ([M)) provides the bound 

Dw(t'm' + l,^'m') < A™ Dw(:/i,l'o)- 

We assume that a > has been chosen large enough so that A < 1 in ([55]) . By estimate ([5^ . it 
follows that 

+ - KA^)\ < K{\\E\\i, ||S||P/')Dw(l^rn' + l,I^rn') ^ 

at every S G 5. So I'm'-i-i converges pointwise to the same limit as I'm', i.e., to i>oo. This allows to 
pass to the limit also on the left-hand side in ((39|) . 

Finally, we show that i/oo is the only solution to ([25)1 . If ly'^ is another solution, i.e., ly'^ = 
T[u'^,u'^], then 

D>v(i^oo,i'L) < ADw(j^oo,i^D- 
Since A < 1, this means that Dw(i'oo, i^^) — 0, and so i^oo = t/'^. □ 

3.2. Symmetries and moments of the stationary state. Having proven the existence of 
stationary state fioo, we now turn to study some of its properties. Our first result is that spatial 
symmetries in the distribution of the coefficient matrices L and R induce symmetries of fioo- 

Proposition 1. Let F C SO(R'') be a matrix subgroup that leaves {L,R) invariant under conju- 
gation, i.e., for any O S F, 

{QLQ-\QRe-'^) = {L, R). (41) 

Then the fixed point ^^o is T-invariant, i.e., pLao o = /^oo for every O G F. In particular, fioo 
is always point symmetric, i.e., /ioo(^) = Moo(^^) for all measurable sets V C K'^. 

Proof. Let X, X' and X" be independent random vectors with distribution /ioo, and let O G F. 
Then, by definition of /ioo as fixed point of Q^, by linearity of O, and by the symmetry relation 

ex = e{LX' + RX") = QLQ-^QX' + QRQ-^QX" ^ LOX' + ROX". 

This means that also Law(8X) — /ioo°0~^ defines a fixed point of Q+. Clearly, /ioooO^^ G Cq{M.'^) 
since /Zoo € Cq(K.'*), and \ldv\ — \v\ because O G SO(M'^). But since Q+'s fixed point is unique in 



Cq, the first claim follows. 

The second claim is a trivial consequence of the fact that one can always choose F — {1,-1}. 
Indeed, one even has {010-^, SRe-^) = {L, R) for = ±1. □ 

The following two propositions are concerned with the finiteness of moments of /ioo- We em- 



phasize that we do no longer assume that p < 3, but we allow for general p > 2 in Assumption 

m 

Proposition 2. The pth absolute moment of ^oo is finite. 
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Proof. This requires just a minor modification of tlie proof of Lemma [H Namely, we estimate the 
remainder term in (1371) in a slightly different way, using 

„/o . 1-2/p 

E[{tT S"r^^-'] < E[(tr5")^/'] ^ < dP^^-'uJ sup 



This means that (1551) is replaced by a different estimate, 

sup — i^^^ — j—f < F sup — ^ , sup 



with the nonlinear function 

The sequence of iterates i^m in the proof of Theorem [2] thus satisfies the following variant of (|55)) . 



<S7max(||S'o| 



\ 1 i^p 

which shows that the p/2-th moment of Vm is m-uniformly bounded. Thus, the p/2-th moment of 
is finite. Now, by definition of fioo in (|26p as a scale mixture, it follows that 

klVoo(d^;)=/ E[|«;^5«;r/2]— -^du;<E[||5||P/2] / \w\P-——dw, 

which is finite. □ 
Proposition 3. Suppose that, for some s > p > 2, 

K* := inf E[|Le|'*] + inf E[|i?e|''] > 1. (42) 

|et = l |e| = l 

Then the sth absolute moment of fioo is infinite. 

Moreover, if {L, R) is invariant under conjugation of the full SO(M'^), then 

< =E[|Leon+E[|i?eor], (43) 

where Bq £ is an arbitrarily chosen unit vector. 

Proof. Let X, X' and X" be a independent random vectors with distribution /ioo- By definition 
of fj,oo as fixed point of , and by the point symmetry stated in Corollary [T] above, 

X' = LX' + RX" = LX' - RX". 

Now, to prove divergence of the sth moment, assume that, on the contrary, E[|X|''] < 00. Then 
the following calculation would be legitimate: 

E[\X\'] = i{E[|LX' + RX"\'] + E[\LX' - RX"\'] } 

i{(|LXf + |i?X'f + 2{X'fL^RX"y^^ + (|LX'|2 + \RX"\^ - 2{X'fL^RX"Y^^} 



E 



> E[(|LXf + \RX'' 

The last inequality follows directly from the convexity of z i-^ z"/^. Observe further that since 
s > 2, one has (a2 + &2)s/2 > a^-\-b^ for arbitrary non-negative real numbers a and b. Consequently, 
with random unit vectors e' = X'/\X'\ and e" X"/\X"\, 

E[\X\'] >E[\LX'\'' + \RX"\'] 

= E\\X'\'E[\Le'\'\iX',X")] + \X"\'' E[\Re"\'\{X' , X")] 

> <n\m- 

This is a contradiction, showing that our assumption E[|X|-'] < 00 cannot be true. 
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To prove formula (|43)) . first observe that for any two unit vectors ei,e2 G S'*^^, there exists 
some Q G SO(R'*) with 62 = Qei. Thus, since {L,R) are invariant under conjugation by Q, and 
Q preserves the norm, it fohows that 

E[|Lei|^] = E[\eLQ-'e2n = E[\Le2n 

and hkewisc for R. □ 

3.3. Regularity of the stationary state. Below, we show that — under rather natural condi- 
tions — the stationary state fioo possesses a density of a certain Sobolev-regularity. 

Lemma 9. Under Assumptions 1-3, let he defined by (j25p . Assume that P{L^e = R^e = 0} = 
and V{L'^e 7^ 7^ R^e] > for every unit vector e. Then the function p 1— > sup|g|^]^ ^(ps) 
decreases monotonically to zero. Consequently, the stationary distribution is absolutely continuous. 

The example in Section 16.31 does not meet the hypotheses of this lemma, and, indeed, its 
stationary distribution is not absolutely continuous. 



Proof. Define the auxiliary functions $, Qx : 

eA(e) *(Ae) 



by 



$(e) lim eA(e). 



(44) 



Recall that ^ is a mixture of Gaussians, thus the function Q\ is real, non-negative and non- 
increasing in A S R+ for every fixed e. Therefore, the limit $(e) G [0, 1] exists for every e. As the 
decreasing limit of the family Qx of continuous and bounded functions, $ is upper semi-continuous, 
and therefore attains its supremum, say at e*, i.e. 

$(e,) = e := sup $(e) (45) 

e|=l 



for some e* e S''-^ Observe that limA^oo *(AL-' e*) = $(£-' e^/jL-' e*|) for all (L, i?) with 
L^e* ^ 0; and similarly if i?-^e* ^ 0. Denoting by 1a the indicator function of a set A, by 
dominated convergence one can write 



= lim ^'(Ae* 



A 

= E 
< 



lim E[^'(AL^e,)*(Ai?^e*)] 

A— ^00 



L e* 



■E 



R'e, 
iR^e, 



E 



R^ 

\R^e, 



92 + (1 _ q)e, 

with q := PfL-'^e* 7^ 7^ R^e^,}, which is positive by assumption. Hence, either 
(which we need to exclude). However, 

1 = $(e,) = *(e,) = E[exp(-i ef S'e,)] 



= or 6* = 1 



(46) 



if and only if eJSe^, = a.s. But this means that — e^E,[S]e = e^S^e,, contradicting the 
hypothesis that S* is positive definite in Assumption [3] 

To conclude the proof, define e„ e § as one maximizer for Ox^ for, say A„ — n, i.e., 

6n(e„) = maxe„(e). (47) 

By compactness of §, we may assume w.l.o.g. that e„ — > e' G § as n — >■ cxd. Now, for given 
e > 0, there exists an = A^(e) such that. 0jv(e') < e/2, and there exists a 5 > such that 
|8Ar(e) — 8jv(e')| < e/2 whenever |e — e'| < S. Since 0„ converges monotonically in n, it follows 
that 0„(e) < e for all n > A^ and |e — e'| < 6, and thus also 

e„(e„) < e (48) 

for all n sufficiently large. □ 

The following is our main result on the regularity of /ioo- 
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Theorem 3. Let Assumptions 1-3 hold. Assume further that there exists a 5 > Q such that 

m := sup E[niin(|L'^e|, |i?^e|)~*'] < oo, (49) 

|e| = l 

and that there exists ana>S such that 

M := sup E[max(|L'^e|, |i?^e|)""] < oo. (50) 

|e| = l 

Then, for every a € (0, a), there exists a finite constant C{a) such that 

\^m<c{am-'^ (51) 

holds for all ^ 7^ 0. In particular, if a > d/2, then the density of fi oc belongs to the Sobolev space 
H'iW^) for any s < a -d/2. 

Proof. Without loss of generality, we assume m > 1 and M > 1. Let a E (0, a) be fixed. Because 
of condition Lemma IHl applies. Thus, we find a p > such that 

l*(e)l<^ forall|e|>p, 

with M defined in dSO]). Next, choose d > such that m-d^^'^ < 1/{3M). By Holder's inequality, 
and in view of this guarantees that 

E[| max(|L'^e|, |i?^e|)-''l„in(|ire|jflTe|)<^] 

< (E[|max(|i^e|, |i?'^e|)-^] )"/>{min(|L^e|, |i?^^e|) < i?}!-"^/" < Mmd^/^ < - 

3 

for every unit vector e. Finally, let C(a) := 2>{p/-dY. 

Now, fix a vector ^0 7^ in W^, and let (Lk,Rk) for k — 1,2,... be i.i.d. copies of {L,R). 
Accordingly, the cr-algcbra Bk is the one generated by {Li,Ri) with 1 < £ < fc. Define random 
elements £,k,'nk G K"^ inductively as follows. Given S,k-i with fc > 1, consider the two random 
vectors v = Lj^^k-i and w = Rj^^k-i, and define 

& = 1|„|>|„,|W + 1|„|<|„,|W, rik = 1|„|<|^|U + l|„|>|^|ii;. 

Thus, ^k is always the "larger one" of v,w. Further, introduce the random variable Xk by 

Xfc = l/(3A/) + l|„,|<^l^,_^|. 

Two properties of Xk will be important in the following: 

K[xk\Bk-i] < ^ and E[|ar"Xfe|Sfe_i] < ^l^-ir^ (53) 
The first part of ([55)) follows directly from the definitions of Xk and since 

E[xk\Bk-i] =^+H\Vk\m-i\<^\Bk-i) 

1 12 
< — + sup P(min(|Lre|,|i??'e|) < t?) < — + m-d^ < —. 

The second part of ([5^ follows with the help of ([5^ . since 

E[\^krxk\Bk-i] < E[(i6i/ia-ii)""xfc|^fe-i]ia-ir" 

< supE[max(|L^e|,|i?^e|)-''(l/(3Af) + l„,i„(|LT,|,|flT,|)<^)]|a._i|-'' 

|e| = l 

As an intermediate step, we are going to verify that 

\^{^k-i)\ < E[\^{^k)Xk\Bk-i] + ^la-il- (54) 
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holds almost surely. Observe that inequality (l54l) is trivially true if < pj'd by our definition 

of C(a), thus we may assume |^fc-i| > p/i?. Since ^I^ satisfies 

|vi'(a.-i)i = |E[vi/(ii^a-i)*(i?ra-i)|Sfe-i]| <E[|vi'(a)ii*fe)i|Sfc-i] 

<E[|M'(Cfe)|((l/3M)l <E[|vl/(efc)xfc|6fc-i]. 
This shows From here, we derive by induction that 

|vI/(eo)| <E[|vI/(^,)|x,x^_i...Xi] +C(a)(l-(2/3/)|eor^ (55) 

Indeed, for ^ = 1, the claim ([55)) is identical to (j54p with fc = 1. Now assume that ([55]) holds for 
some £ > 1. Inserting ((54|) with fc = £ + 1 into (j55p gives 



|*(^o)| < E[|E[|vl/(e,+i)|x£+i|S,] + ^\^,\-\x,Xi-i • • • Xi] + C(a)(l - {2rS)%or 

< E[|4'(C,+i)|x^+iX£ • • • Xi] + ^IE[ie£r"x^X^-i • • • Xi] + C{a){l - (2/3)0 leol"". 



Apply the second estimate from to the second expectation above 

2 
3 



E[ied""X£X£-i ■ ■ • Xi] = E[E[|ed-'^x^ |^£-i]x<'-i • • • Xil < ^E[|^,_i|-''x£-i • • • Xi] 



<---<{2/3)%r. 

This verifies (j55p with ^ + 1 in place of £. To finish the proof, observe that (l55|) for every £ > 1 
implies ([5T|). since, by (|53)) . 



0<E[|«'(e^)|xm-i---Xi] <E[E[x^|e,_i]x^-i---Xi] < ^E[x£-i---xi] 

< • • • < (2/(3Af))', 

which tends to zero as £ tends to infinity. 

The statement about Sobolev regularity of the density function follows since if the characteristic 
function 4" satisfies (j5ip . then |CI*^(0 belongs to L^(R'') for every s < a—d/2. In turn, the inverse 
Fourier transform of ^ is a density function in the space H^(M.'^). □ 

4. Convergence to equilibrium 

Having established the existence of a stationary distribution for ([I}, we now turn to study 
transient solutions and prove their convergence towards equilibrium. 

Transient solutions fj,(t) to the initial value problem ^ are easily obtained, for an arbitrary 
initial probability measure /xq on R'^, by the classical Wild construction [53]: define inductively 
probability measures fin by 



1 " 



or, equivalently — recalling — in Fourier space, 

An+i(e) — Y5]E[A,(i^e)A«-,(i?^0]- (56) 

j=o 

Then the infinite convex combination 

oo 

defines a weak transient solution to ((T|): the curve t i->- fi{t) is weakly continuous with respect to 
t > 0, and satisfies 

^ [ ^{v)fi{t;dv)+ f <j,{v)fi{t-dv) ^ f (t>{v')Q+Wmv') (58) 

It JRrf JK,<1 JRd 

for arbitrary test functions (p G Cb(E'*), and at every t > 0. 
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4.1. Qualitative results by probabilistic methods. In this section, we apply probabilistic 
methods related to the central limit theorem to prove the following qualitative result on conver- 
gence to equilibrium. 

Theorem 4. Assume that the initial condition has unit temperature, 0[/io] = 1- Then the distri- 
butions fin converge weakly towards the stationary state fioo obtained in Section\^ above. Conse- 
quently, the associated transient solution fi{t) converges weakly to fi^o as t ^ oo. 

Our key tool is a probabilistic representation of the Wild sum ([57]) . following [3J [TH]. On a 
sufficiently large probability space (f2, J^, P), let the following be given: 

• a sequence (X„)„gN of i.i.d. random vectors with distribution /iq; 

• a sequence ((L„, i?„))^^^j^ of i.i.d. random matrices, distributed as {L,R); 

• a sequence (/„)„gN of independent integer random variables, each /„ being uniformly 
distributed on the indices {1, 2, . . . , n}; 

We assume further that (/n)neN, {Lm Rn)neN and (X„)„gN are stochastically independent. Define 
a random array of weights [/3j^„ : j — 1, . . . ,n + l]„>o recursively: Let /3i_o ■— 1, (/3i,i,/32,i) ■— 
(Li, and for any n > 2, 

Wl.n, • ■ • , /3n+l,n) := (/3l,n-l, ■ • ■ , /?/„ -l,n-l , /3/„ .n-l^n , Pln,n-lRm /3/„ + l,n-l , ■ • ■ , Pn,n-l)- (59) 

Finally, setting 

n+l 

:=^/3,-„X,, (60) 
i=i 

one obtains the following alternative representation of the measures jin from (|57p . 

Proposition 4. For every n>l, one has 

A„(e) =IE[exp(zC^Ty„)]. (61) 

Proof. The proof is analogous to the proof of Proposition 1 in 3^ for the scalar case, with the 
only difference that here one needs to be careful about the order of multiplication of the (non- 
commutative) matrices and Rk . We briefly repeat the argument as we shall refer to parts of it 
later. 

We prove by induction on n > 1. For n = equation is clearly true. For ?i > 1, let 
J„ be the random index such that 

n ^ / -^i/Sjn for every j < J„, 

P^''' \ i?i/3;„ for every j = J„ 1, . . . , n + 1, 

with suitable random matrices /3j„ and /3J„- It can be shown — see Proposition 1 of [3] — that 

(i) Jn is uniformly distributed on {1, . . . , n}, and 

(ii) (/3| „)j=i,...,j„ and (/3J„)j=j„+i,...,n-i-i are conditionally independent given J„ with condi- 
tional distribution equal to the distribution of (/3j,j„-i)j=i,...,j„ and (/3j_„_j„)j=i^...^„+i_j„ , 
respectively. 

This shows that 

J„ n+l 

J=l ]=J^+l 

for two sequences (Wq, . . . , W^) and {Wq , . . . , W^) of i.i.d. copies of (Wo, . . . , Wn) independent 
from Jn and {L.R). Hence, using the induction hypothesis, 

^ n — 1 

= E[A,7„-l(L^0A«-J„(i?^0] = —-Y.nl^^{L^O^^n-^{R^O] - A« (C) • □ 

n + l — ' 

As a first consequence of the probabilistic representation we prove 
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Lemma 10. Under Assumption 1, and if ij-q is centered, then fi{t) is a centered probability measure 
for allt > 0. If additionally /io has finite temperature ©[/io] < oo, then this is conserved, 0[/i(i)] = 
0[/io]. // instead 0[/io] — +oo, then also &[^{t)] = +oo for all t > 0. 

Proof. It is clear that it suffices to prove the statement for First of all note that (O yields 
-^[ll-^lli + ll^lli] < hence i?[||_L||2 + II-RII2] < +00. As a consequence, we shall show that 

Sj=i /^jnlh] < +00. It clearly suffices to show that E[||/3j„||2] < +00 for every j = 1, . . . ,n + l. 
Now each /3j„ is a random product a random number, say Sj, of matrices L and R. Since (L^, Rk)k 
are independent and identically distributed and, for each fc, Law(Lfe) — Law(i?fe), using the 
subadditivity of the norm and the fact that Sj < n a.s., we have 

nmnh] < nil WLtU < nil max(l, WlkM < (E[max(l, IIL1II2)])" < +CX). 
fc=i fc=i 

Propagation of centering follows now from the linearity of the collision rules ([3]) . Indeed 

n+1 n+1 

E[I^„] = pj^Xj] = E[^ /3,n]E[Xi] - 0, 
j=i j=i 

since {Xj)j>i and (/3jn)j"=i,...,n+i are independent and E[Xi] = 0. As for the temperature, 

n+1 n+1 

Clearly, since the Xj are i.i.d. centered random variables, and are also independent of the /3's, 

n+1 n+1 n+1 

^ E[Xj/3j„ A„X,] =0, and ^ E[Xj/3j„/?,„X,] = E[Xf E[^ PfnPjnjX^] . 
■'.J=.' J=i J=i 

It remains to prove that 

n+1 

nJ2PfnPjn]^l- (62) 
J = l 

We shall prove (EH) by induction. Set On ■= Y.]=l ^JnPjn- We have E[6I„] = 1 for n = and 
n = 1. For n > 2, with the same notation of the proof of Proposition 21 we can write 

Jr^ n + 1 

E[0„]^E[^(Li/3W)^ii/?S+ E (^i/^S)'^^i/5]n] 

J=l J=J„+1 

,/„ n+1 

i=i j=J„+i 

Since E[i?fi?i] — E.[LfLi] = ^1 and and /3j2)/^jn'' ^'"^ stochastically independent, recall- 

ing also (i) and (ii) of the proof of Proposition 01 we can write 

Jn n+1 n 

i=i j=J„+i fc=i 

Hence we can conclude by the induction hypothesis. □ 

Theorem [H will be proven by studying the distributional convergence of the random variables 
Wn with methods related to the central limit theorem. For definiteness, denote by B the cr-field 
generated by the weights /3n,j, for n > with 1 < j < n + 1, respectively. Accordingly, introduce 
the conditional covariance matrices 

n+1 

Sn := E[WnW^\B] = E l3n,jSof3l, and := Law(5„), (63) 
j=i 
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SO that in particular 5*0 = Cov(/io) a.s., and vq := Ss^. In a series of lemmas, we show the weak 
convergence of the i>n towards v,^ obtained in Theorem [51 First, we derive an a priori estimate 
on the weights. 



Lemma 11. Under Assumptions^^ and\^ there exists a constant C such that 

n+l 

M(") := ^ sup E[|/3j„er] < Cn^^-' 



(64) 



for every n > 1 . In particular, Mp — )■ as n oo . 
Proof. We first show a similar estimate for a related quantity, namely 

Here F denotes Euler's function. We proceed inductively with respect to n. By construction of 
the weights, we have 



(65) 



j)#0 



r//0 W(7?) 



a;(r/) 



(66) 



Since (L„,_R„) and Pk..n-i are independent, it follows that 



E 



uj{r]) 



Br, 



< E 



uj{r)) 



sup E 



uj{r]') 



< KpE 



where Assumption [5] and the symmetry {L,R) = {R,L) have been used in the last step. Thus, in 
we have 



supE 

>)#0 



sup 1 



a;(r7') 



< Kp sup E 

»,#0 



r'^(^L-i'7) 



a;(ry) 



which leads to the recursive relation 



and implies (|55|) . The original claim ([M)) now follows immediately, since 

Lemma 12. Under Assumptions 1-3, i/ 0[/io] = 1; o^c' with he defined in Theorem\^ it 
follows that 



r(") 



Dw( ) < - / , Dw('^'i-i, i^o 



(67) 



where Dyy is t/ie metric introduced in (j29p . to^/i a Zan7e enough such that X < 1 in psp . Conse- 
quently, the Vn converge weakly to Voo- 



Proof. Since is a solution to (1281) . we can write 

1 

n — ' 

i=l 

Recall the definition of t/„ in (|63p . We prove that for every n > 1, 

1 " 

77 ^ / 



(68) 



(69) 
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Indeed, with /3| „, and J being defined as in the proof of Proposition |4l we have 

J n+l 

Now follows since the expressions in the round brackets are independent of each other and 
of (ii,i?i), and are distributed like Sj-i and Sn^j, respectively. Using the representations 
(|68 )l fc ((69| and the convexity of Dyy, see ((3T|) . we obtain 



1 " 

Dw( 

77 Z ^ 



Now ([ST]) follows from p4l) . with A < 1 defined as in ([55)1 . To conclude the proof, we recall that 
if a sequence (a„)„ of non- negative numbers satisfies the inequalities 

A 



77 Z / 



for every n > 1, with some fixed constant A G (0, 1), then 

1 - A 

a„<ai[](l _) <Cai7i-(i-^), 

which tends to zero as n — > oo. Inequality (j67p thus implies D>v(j/„,i^oo) — ^ 0, proving weak 
convergence of to i/qo- D 

Weak convergence of the Vn is not quite enough to apply the methods of the central limit 
theorem. We also need a uniform control on the second momenta, which is related to the Lindeberg 
condition. 

Lemma 13. Under Assumption 1-3, and i/0[/io] = 1, it follows that for every e > and every 



.n+l 



hm E\j2\ef^MH{\ePjnX,\ > e} 



0. 



Proof. Let ^ e M'' be fixed. For notational simplicity set 



(70) 



(71) 



Since E[|Xj p] = d by hypothesis, for every 77 > there is a i^T = if (77) such that E[|Xj pI{|Xj | > 
K}] < 77. At this stage observe that, by ([5^ . 



_ n+l 



n+l 



J = l 

Moreover, the rule (fT7|) implies 
and so we obtain 

n+l 

j=l 



=E[^tr(/3j„/3,„) 



= trl = d. 



i=i 



\X,n\^ < mn\\l\^\^\X,\' and |X,„|' < \PU\^\X,\', 



n+l 

[Y,\X,nfl{\X,,,\>e}] <E[Y,mn\\l\\eE[\X,\H{\X,\ > K} 



(72) 



(73) 



j=i 



_ n+l 



E 



[j2\Xjnm\Xjn\>e,\X,\<K} 



i=i 



n+l 



< 



dv\e+K'E[T.\pfneH\PU\>j^} 



J=l 
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Let q p/2 and q' := p/{p — 2) be its conjugate exponent. Then, by Holder's and by Markov's 
inequality, 



E 



\pun{m\ > ^}] < {mnmr {^[\pu\ > 



Hence 



. n+l 



E 



n+l 



r;d|eP+X2(f)''|CrM("). 



With dM]), we finally get 



n+l 



limsupEf^ |X,„|2l{|X,„| > e}j < 7yd|Cp 



Since 77 is arbitrary, this shows (1701) . 

This concludes our preparations for the proof of Theorem SI 



□ 



Proof of Theorem^ Let be a random variable with law Urx,. Let {Yj)j be a sequence of i.i.d. 
random vectors with Gaussian distribution of zero mean and covariance matrix Sq. Observe that 
Sn = GoY\Y^^t:i PjnYj\B]. Since a linear combination of Gaussian random variables is a Gaussian 
random variable, 



E 



n+l 



expK^^(/?,„r,)} 



= E 



= E 



E 



n+l 

cxpK^^(/3,„y,)} 
1 



B 



1 



exp{--(^^5„0} exp{--(C^5o,0} 



Here we used the distributional convergence of Sn to 6*00 by Lemma [121 Thus, Theorem |4] follows 
if we prove that 

n+l n+l 



E 



^M^e Y.^PjnXj)}\ - E[exp{je^ Y.^PjuYj)) 
i=i i=i 



^ 0. 



(74) 



Define Xjn as in (j7ip above, and accordingly Y,„ :— {f3jnYj), as well as the conditional second 
momenta in direction of ^, 

:= E[Xl\B] = E[Y^JB] = eiPjnSoPfjC 
Also, we introduce the characteristic functions 

(j)jn ■■= E[exp{iXjn}\B], ijjn ■■= E[exp{iYj„}|6]. 
Clearly, ((74)) is equivalent to 



_ n+l 



n+l 



E 



0. 



Since for arbitrary complex numbers zi, . . . , z^i and z[,...,z!^ of modulus less or equal to one, 



j=i j=i i=i 



(75) 
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we can write 



< E 



-l + -<7 



E 



-l + -a 



(76) 



By formula (jl02p from the appendix, any centered random variable Z satisfies 



L[e'^] - 1 + ^E[Z^] 



< E 



(77) 



Observing that 

E[Yj,,\B]^E[Xj„\B]^0, 
we can use ((7S)l - (|77|) and the definition of (pjn and ipjn to conclude 

n+l n+1 

|E[ Y[ (j)jn - n V^J"] I < M^) + Bn{e) 
j j 

for arbitrary e > 0, where 

^«(^) ■■^\Y.^mn?l{\X,r.\ < e}] +^E[|X,„pI{|X,„| > e}], 
j j 
and Bn{e) is the respective expression for Yj„. In view of (j73p . we can estimate 

An < J|CpE|Xip^E[||/3,„||2] +^E[|X,„pi{|x,„| > e}]. 

j j 

The first term on the right-hand side can be made arbitrarily small for an appropriate choice of 
e > 0; recall that E[||/3j„||2] = d hy ((72)) . And the second term converges to zero for n — >■ cx) 
by ([70]) . independently of e. The analogous argument applies to i3„(e). □ 

4.2. Quantitative results in Fourier metrics. Under a very weak additional hypothesis on 
the initial condition, the previously obtained qualitative convergence result can even be quantified 
on the level of the Fourier transform. 

Theorem 5. Let Assumptions 1-3 hold, and let he positive definite. If the initial condition jj^o 
belongs to a class C2+a(^'^) with some a > 0, i.e., if its absolute moment of order 2 + a is finite, 
then there exists a constant A € (0, 1) that depends on /iq only through a , such that 



(78) 



|A(^;C) - Aoo(C)| <C(/io)max( 
holds for allt>0 and all ^ G M'*. 

The appropriate Fourier distance in the current context is given by 

dUl^i,fi2) = sup {uj{0''\fii{0 - A2(0 + ^eiCovifii) - Cov(m2))C|) 

+ a|i Cov(^i) - Cov(/X2)||oc., 

with a parameter a > that needs to be determined. This definition is similar to the one for Dw 
given in (j^ . 

Lemma 14. The distance defines a metric on Cp(M'') with the following convexity property: 
// /ii — s^'i + (1 — s)/ii and /i2 — s/i2 + (1 — s)/i2 are two convex combinations of measures in 
CpiR'^), with s e [0, 1], then 

dUl^i,fi2) < sd^{^l[,^I'^) + (1 - s)d„(Ai'i',A*2)- (79) 

Moreover, d^^ provides a pointwise control on the Fourier transforms, 

|Ai(e)-/i2(0| <i^niax(|C|Mer)d^(/.i,/X2), (80) 



24 



FEDERICO BASSETTI AND DANIEL MATTHES 



with a constant K that only depends on a and uj. Finally, convergence dt^{^k, fJ-*) ^ for a 
sequence {^k) in CpiW'') implies weak convergence of fXk to /i* and convergence of the covariance 
matrices. 

We omit the proof, which is very similar to (and easier than) the proof of Lemma [5] 

Lemma 15. Let Assumptions 1-3 hold. Let /ii, /i2 G Cp{W^). Then 

d^(Q+[/ii],Q+M) < Ad^(/ii,/i2) (81) 

with the same constant A defined in (|35p . 

The proof combines the elements from the proofs of Lemma [5] and Lemma [T] We leave it to 
the reader to adapt the technical details accordingly. 

Proof of Theorem\^ Without loss of generality, we may assume that a < 1 and p' := 2 + cr < p. 
By Lemma [1] Assumption [2] is satisfied for p' with a suitable weight function lo' . Now apply 
Lemma 1151 with lo' instead of oj. Choosing a > sufficiently large in the respective metric d(^/, 
one obtains 

d^'((3+[/ii],Q+[^2]) < Adi^'(/ii,/i2) with A := max (k^/, k + 2a"^d^'/^Kp/) < 1. 
Combining the convexity property (j79p with the fact that, for any s > 0, 

(Q+[/^(s)],Q+[Aioo]) <Ad^'iKs) 

it follows that 

dtj' (/i(i), /ioo) < e *dtj/(/io, 

(Q+[m(-5)] 1 Moo 

)ds 

Jo 

5 Moo 

)+A / e-(*-^M<.,(/i(s),/ioo)ds. 
Jo 

Gronwall's lemma now yields 

d„'(Ai(i),Aioo) < d^'{^io,^Ioo)e-^^~^^'. (82) 
From dMl), the claim dTS]) follows with the help of dHHl) □ 

5. Divergence for data of infinite energy 

In this section, we prove that the solution to the initial value problem ^ cannot converge to 
a stationary state if the temperature of the initial condition is infinite. For z e M'' and r > 0, let 

< r} be the closed ball of radius r around z. 

Theorem 6. Let Assumptions 1-3 hold, and suppose that the stationary solution ^irx, of unit 
temperature has no concentration in zero. If fi{t) is the solution to ^ with 0[/io] — +oo, then 
/.t(t)[Br(0)] converges to zero as t ^ oo, for every r > 0. 

In the following, let Z > be a parameter, which is large enough such that 

ai :=/io[B,(0)] > 0. 

Accordingly, we decompose /io = aj/i^ + (1 — a;)/Io, where the probability measures /i^ and JIq are 
given by 

Mo[-nB;(o)] _ Mo[-nB;(o)^] 

n[-\ := and /iQ ' ■= 1 , 

ai 1 - a; 

and introduce expectation and variance of the restriction. 



mi := / vn{dv) and o", := / {v - mi) nJdv). 

Finally, recall that 5* — /too is the characteristic function of the stationary solution with covariance 
— see ipS]) — and define 

H{u):= [ e-'^^{^u)d^^^ [ e-^^i^d^. 
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Observe that H{u) — > as m — > cx) if and only if fioo has no concentration in 0. The key is to 
estabUsh the fohowing time-dependent upper bound on ^{t)[Br{z)]. 



Proposition 5. Under Assumptions 1-3, and if r > and ai > \plr , then 

sup MW[Br(^)] < c[(Z/r)f max(e-(i-^)*,e-i(i-''^)-*) + i7(V^c7i/ryd)' 



(83) 



where A e (0, 1) is the constant defined in Theorem\S\ for 2 + a — p. 



Proof. For the sake of notational simphcity write a for ai and m for to/. Define the array 
{/3j,n}n>o,i<j<n+i of random weights as in ([5^ . In addition, consider a sequence of indepen- 
dent and identically distributed random elements 

where aj,X^j, Xj are independent of each other, and independent of the Pj^n- The distributions 
of ttj, 2Lj and Xj, respectively, are given by 

P{aj = 1} = 1 — P{aj = 0} ~ a, I^q — LawX^, JIq ~ Law Xj. 

It is easily seen that Xj :— aj2Lj + (1 — <^j)^j has law fiQ, that \2Lj\ < I < \Xj\ a.s., and that 
TO = K[Xj] and af = E,[\2£j — to^], for all j > 1. In particular, in the Wild representation (1571) of 
the transient solution /i(t), we can write /i„ = Law[^jl^-j'^ /3j^nXj]. 

Further, define the Gaussian density 7(2/) = exp{—\y\'^/2). By Parseval's identity, we obtain for 
any r > and z in R'* that 



(2tt) 2 Jmd 



Since ^{{v — z)/r) > e ^1"^ for all v G B,.(z), it follows further that 



(27r) 



62 



i- supAi(t)[B,(z)] < (27r)2 / 7 



V — z 



Kt;dv)< / j{om;^/rm=--nt). (84) 



Denote by (?!)(^) the characteristic function of X.i — m and by $ the characteristic function of Xi. 
For the sake of simplicity set r„ = e~*(l — e~*)". The Wild sum representation implies 



_ ri+l 



ri+1 



n+1 



n>0 



J = l 



n>0 



For further estimation, introduce the semi-definite matrix as the square root of S/ = Cov(X^). 
Clearly, tr(I]/) = af. An application of the triangle inequality yields 

I A(^; I < E ^" (^1." + + ^3,n) 



rt>0 



with 



n+1 



Ai,„(0 := e[| n ^(a,/3j:,0 - n7(".Vs;/3f„0 



n+1 



^3,n(0 := n7(V^/3j„0- 
Since the last term is real and positive for all ^ G M'*, no absolute value needed. 
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Estimate for Ai_„. Another application of the triangle inequality and of ([75]) gives 



n+l 



n+l 



n+l 



The differences inside the absolute values are now estimated using (I102p from the appendix. Since 
7 is the characteristic function of a Gaussian random vector Y (with zero mean und unitary 
covariance matrix), this leads to 

n+l 

Recalling estimate (j64p on the sums of /3's, we arrive at 



A, 



l.^{l) <Ciiaf + mrn 



Estimate for A2,n- Since |e 



where r„ := 



r n!r(Kp) 
e « I < l^^ — y^l and since p/2 > 1, we have 



< E 



An application of the von Bahr-Esseen inequality (I103p reveals that 

\A2A0\^ < ^nai-af^J2^[\ei3jn^il3U\^ < ^^Ejai - a|i E sup E|/3jeP 

and finally yields 



\2/p 



Estimate for A^^n- Define n*{t) for t > by 



_ n+l 



n>0 j = l V ' 

which is the solution of the Boltzmann equation with initial datum 



■(!)({) ^=E-»'^[nHi/^Cf 



In particular, each has mean zero and unit temperature. Thus, we can apply the quantitative 
estimate ([75]) for convergence of ^*{t) to /ioo, obtaining 



_ n+l 



_ n+l 



Y.rn^[\{l{V^.PUlr)] < T.^nE[l[ifjpfn^C 

n>0 j=0 n>0 3 = 1 ' 



< Cd, 



,(^*(i),/Xoo)max 



2 2 



'aa 



\fd ' r\/d 



< C'e-"-«l' max (i, ^) ma.x(|ir, IJR + 
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Combining the estimates on Ai_„ to A^ ^i yields the following bound on I{t) defined in (j84p 



n>0 

Using the fact that 



g -(l-Kp)t 



n=0 

and that — since r„ > with X]^o ~ since z i— >■ z^/p is concave 

E r„r2/f < ^ r„r„ = exp ( - - (1 - , 

11=0 ^11=0 ' ^ 

we conclude that 



sup 



Since = E[(X — m)^] < 2P, and under the assumption that ct/ > r-\/2, we eventually arrive at 
formula □ 

Proof of Theorem\^ It remains to be proven that estimate (j83p implies Theorem [6l To this end, 
observe that for every e > 0, there is / such that for every I > I 

.2 ^ 1 - e f L ,2, 



See, e.g., [7]. Hence, if 0[/io] — +oo, then also ^Jaloi — > +oo as I — > +cx). Consequently, the 
second term in the sum (I83p can be made arbitrarily small by choosing I large enough. And the 
first term in the sum converges, for fixed to zero as t — cx3. □ 

6. Examples 

6.1. Example: Maxwell molecules. In this section, we discuss one example for ((U, which 
provides a model for a particular case of inelastically colliding Maxwellian molecules [S] ; see [TU] 
for a recent review. The linear collision rules are given by 

v' ~ V ^ an • (u» — u)n, v'^= — an • (u* — u)n. (85) 

These collision rules are brought in the form ^ by setting 

L = ann^ and i? = 1 — ann"^. (86) 

Physically, the unit vector n points along the contact line of the molecules upon collision. In 
accordance with our basic model assumptions, n is supposed to have a fixed distribution on the 
sphere that does not change with v and w*. For simplicity, we shall consider a uniform 

distribution, but situations with spatial anisotropics could be discussed along the similar lines. 



Remark 2. The rules (j85p are given in the so-called ui -representation, which can easily be 
rephrased in terms of the (standard) a -representation, see e.g. jSlj, 

V = + = — —<r. (87) 

For Maxwell molecules, one typically assumes an easy distribution for a on W^^^ , referred to as 
cross section, which leads, after the corresponding change of variables, to a complicated distribution 
for n. In turn, the uniform distribution for n assumed here gives rise to the cross section 
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The non- negative random variable a in (1851) models the inelasticity of the collision. Our model 
is thus similar — but not identical — to the one for inelastic Maxwell molecules with a background 
heat bath proposed in [TT]. For a = 1, one obtains a particular model for elastic Maxwell molecules, 
with the untypical cross section (|88p . The resulting rules (still for a = 1) are very strict in the 
sense that they conserve the particle momenta and kinetic energies in each individual collision, as 
becomes evident from 

L^L + R^R^l a.s. (89) 

Indeed, if one chooses the strict correlation {L,R) = {L^,,R^,) in Q, then 

V* + vi = {L + R)v + {L + R)v^ =v + v*, 

and, since L'^L = L, R^ R = R and R = R^ L = 0, 

|t,'|2 + |^,;|2 ^ v^L^Lv + v^R^Rv, + v^L^Rv, + v^L^Rv + |<p 

= v'^ {L + R)v + vj {L + R)v, = |fp + \v,\^. 

The stationary state ^oo is given by the Maxwell distribution d//oo/dw = (27r)~''/2 gxp(— 
This follows, for instance, since (15^ implies that S = 1 a.s. is a solution to ([25)) . In the general 
"inelastic" case, where a is not identically one, we require in addition that a is a positive random 
variable that is independent of n and satisfies 

E[a(l-a)]=0. (90) 

For simplicity, we also assume that a is "concentrated near one" in the sense that 

E[a^{l- af] < cE[a^] (91) 

for a constant c > determined below. 

In preparation of the following calculations, introduce the random variables Vi to Vd by 

V, = ef n, (92) 

where n is uniformly distributed on the sphere S^"-'^ C W^. Notice that Vi + ■ ■ ■ + — 1 almost 
surely since n is a unit vector. Thus, 1 = E[Vi + • • • + V^] = dE[Vi] since the Vi are identically 
distributed, and so E[V]2] = l/d. For later reference, we also remark that 

1 = E[{V^ + ■■■ + V^f] = dMlV^] + d{d - l)E[V^V^] (93) 

since each product ViVj with i ^ j has the same distribution as V1V2, from where it follows that 
E[Vl] < 1/d. 

We are now going to verify Assumptions 1-3 for the model at hand. For notational simplicity, 
we shall directly work with the unsymmetrized matrices (|86p . 

6.1.1. Assumption[Ii Since nn-^ defines an orthogonal projector for every n € S''"^, and recalling 
that n and a are independent, it follows that 

E[i^i + R^R] = 1 - 2E[a(l - a)]E[nn^]. 

Condition ((90t yields the desired result. 

6.1.2. Assumption[^ We are going to prove (fTU|) for the weight function a;(^) = |^|'*. For an 
arbitrary unit vector e, we have 

E[|L^e|'* + \R^e\'^] = E[a*]E[F];*] + 1 - 2E[a(2 - a)]E[V^] + ^[0^(2 - af]E[V^] 
= 1 - 2E[a2]E[T/j2] _,_ 2E[a'^ - 2a^ + 2a^]E[V^] 
= 1 - 2E[a^]{E[V^] - E[y/]) + 2E[a^{l - a)^]E[V^]. 
Assuming that the constant c in ([9T|) satisfies c < E[V^/] < E[Vi] — E[Vi], inequality (fTO)) follows. 
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6.1.3. Assumption\^ We already know that = 1. Let a symmetric matrix M G 5 be given, 
with trM = 0. Since the distribution of (L, R) is rotationally symmetric, we may assume that M 
is diagonal with respect to the standard basis, 

d d 

M = ^^lkekel, with ^fik = 0. (94) 

fc=l k=l 

Now, let some unit vector e g S''"^ be given, and introduce the numbers \i = eje, satisfying 
+ • • • + = 1. We obtain, recalling that a and n are independent, 

d d 
i.j^k—1 i,j^k—l 

At this point, observe that 
Thus, 

d d 

e^E[LML^]e ^ Eia"^] ^ X^HkE[Vi^Vi] + E[a^] ^ A|/ifc(E[y*] - E[Vi^Vi]) 

i,k=l k=l 

d d d 

= E[a^]E[VM] E E +E[a2](E[y/] - E[VM]) E 

i=l k=l fe=l 
=0 

For the evaluation of the analogous term with R in place of L, observe that 

ejRek = ef Gfc - aef nn^efc = 5ik - aViVk- 
Using that E[T/iV,] = (5ijE[V\^] and performing similar manipulations as above yields 

d 

e^E[RMR^]e^ E AjAjAifcE[ef i?efce^i?^ej] 

i,j.k—l 
d 

= E ^^^JM^^kSJk - aS,kE[V,Vk] - aS.knV^Vk] + a^E[V,V,V^]) 

d d d 

= E^'^fe-2EH E ^*^feMfcE[F,T4-]+E[a2] E A,A,^feE[T/,F,lf ] 

k—l i^k—1 i,j.k—l 

d 

= {I - 2E[a]E[V^] +E[a^]{E[V^] - E[V^Vi])]^\l^lk■ 

fc=i 

Adding up, and using that E[a^] = E[a] by (|M|) . we have 

d 

e^E[LML^ + RMR^]e = {l - 2E[a] {E[V^] + E[V^Vi] - E[V^]) } E >'kf^k < k max |^fc|, 

fc=i - - 

with the constant 

K = l- 2E[a]{E[V^] + E[V^V^] - E[V^^]). 
Here we have used implicitly that n > 0, which is true since 

E[V^Vi] < E[V^] and E[V^] = l/d < 1/2, 

and since E[q;] < 1, which follows from E[a]'^ < E[q;^] = E[a]. To conclude Assumption[31 it suffices 
to check that indeed k < 1. But this is a trivial consequence of E[Vf^] ~ 1/d and E[V\'*] < 1/d, 
which in turn follows from (1931) . 
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6.1.4. Heavy tails. Above, we have shown that the fourth absolute moment of /ioo is finite if c in 
(Plj) is smaU enough. Next, we argue that unless a = 1 a.s., there exists a s > 4 such that the sth 
absolute moment of /ioo diverges. In other words, the steady distribution /ioo develops fat tails as 
soon as the collision rules are "a bit inelastic" . 

The goal is to show for some (sufficiently large) s > 4. Observe that, to this end, it suffices 
to prove 

E[|i?ein = ¥.[a'\&[V^] > 1. (95) 

Since we assume that not a = 1 a.s., it follows from (|90p that a > 1 with positive probability. 
Consequently, there exists some e* G (0, 1/2) such that p* := P(a > 1 + 2e*) > 0. Observe further 
that q* := P(V"i > 1 - £*) > 0, and so 

E[a-']E[T/i"] > p*q*{l + 2e*)"(l - e*)^ > p*q*{l + e72)^ 

tends to infinity as s — ?► oo. So (|95|) holds for some large s. 

6.1.5. Regularity. We have seen that outside of the elastic case, the stationary distribution ^oo is 
not a Gaussian, but a scale mixture. Depending on the distribution of a, this mixture might or 
might not possess a smooth density. 

In order to check the applicability of Theorem [31 we need to evaluate the criteria (H^ and ((5(1)) . 
By rotational symmetry, it suffices to calculate the expectations for one particular unit vector, 
say e = ei. For (l49t . we find 

E[min(|LV|,|i?V|)-'] <E[a-^](E[|yir^] +¥.[{1 - Vl)-'/^]) 

= j' \z\-'il - z2)('^-3)/2d2 + J\i_ z^)id-3-S)/^dz^ . 

Here ad is the ration of the respective hypersurface measures of S''"^ and S''^^. The sum of the 
last integrals is obviously finite for as long as 5 < 1. Thus (UHl) holds if E[a^*] is finite for some 
(5 > 0. 

On the other hand, since |L^eip + |i?^eip = a.s., it follows that 

E[max(|L^ei|,|i?'^ei|)-^] > E[(aV2)-'^/2] =2^/2E[a-^]. 

In conclusion, the regularity of /ioo is entirely determined by the largest exponent d for which 
E[a~°] is finite. If, for instance a > e > a.s., then a > is arbitrary, and /ioo possesses a density 
in C°°{R^). 

6.2. Random rotations. Let A, B be random matrices in SO(R''), i.e., A and B represent 
random rotations in W^. Further, let a,/? be non-negative random variables with 

E[a2 + /32] = 1 and Kp E[aP + /J^] < 1 (96) 

for some p G (2,3). For simplicity, we assume further that (a,/3) is independent of {A,B). Now 
define 

L = aA, R = /3B. 

This is a relatively straight-forward extension of a random kinetic model on R to R'^. As already 
pointed out in the introduction, examples of this type might violate Assumption [3] if the rotations 
leave certain genuine subspaces of R'' invariant. We claim that our set assumptions is met if the 
following is true for every unit vector e: the measure /ig induced on the sphere S''^^ through A^e, 

fie{B) = P(A^e e B) for aU Borel sets C 

can be decomposed in the form /Xg = ecr + (1 — e)/ie, where a is the normalized surface measure 
on S"*"^, and /ig is some probability measure, with an e > independent of e. 
Since A, B are rotation matrices, A^ A = B^B = 1, and so 

E[L^i + R^R] = E[a^ + /3^] = 1, 
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which shows Assumption [T] Also, with the same p as in (|96|. 

which is Assumption [5] with cj(^) ~ To check Assumption |21 first note that = 1 satisfies 

Next, let a symmetric matrix M g 5 be given. Then, for every e e E>^~-^, 

e^¥.[LML^ + RMR^]e = ¥.[a^ [A^ ef Y.{A^ e) + [B'^ ef i:{B^ e)] 



= (l-e) / w'^Mwfieidw) +e / u;^Mwcr(dw), 

where the assumed decomposition of iie has been employed. The last integral equals to tiM/d, 
and thus vanishes if M is traceless. This is easily verified by decomposing M in the form 
defining Vi, Vrf as in ([Ml) with w in place of n, and observing that the integral equals to a sum 
over /.tfcEfV;?] for fc = 1, d. This means that 

|e^E[LA/L^ + RMR^]e\ < (1 - e) / ||M||ooAe(du;) = (1 - e)||M||oo, 

proving Assumption [3] Thus, our theory provides the existence and uniqueness of a stationary 
distribution ^oo- 

Once this is known, it is now easy to conclude that fioo is — as one would expect — the 
rotationally symmetric extension of the stationary state /x* € V(R) for the corresponding one- 
dimensional model with scalar coefficients {L'.R') — (a,/3). From the results in [3], it follows 
that the one-dimensional equation possesses a unique stationary solution fi^ that is centered and 
of unit second moment; call its characteristic function i/j. We claim that 5*, the characteristic 
function of fioo, satisfies "0(1^1) for all C G H^''- In fact, this follows immediately since 

Ei^iL'^O^iR^O] = mialA'^mWlB^m = E[7/^(a|e|)0(/3|^|)] = m\) = 

6.3. An asymmetric example. The last example equation we consider produces a steady distri- 
bution that is not rotationally symmetric. Thus, even if rotationally symmetric initial conditions 
are used, the results from [6 do not apply, as the symmetry is broken at any t > 0. Specifically, 
we state the example in the plane M^, and the stationary distribution /loo will be concentrated 
on the coordinate cross {{vi,V2) & M.'^\viV2 — 0}. (It is not hard to generalize the idea developed 
below to higher dimensions and to more sophisticated non-symmetric stationary states.) 

Recall that ei, 62 are the canonical basis vectors, and introduce the vectors e+, e_ by e± = 
(ei ± 62) /V^- The random coefficients L and R are defined by 

where the random index i attains the values 1 and 2 with probabilities q and 1 — 5, respectively. 
Since e_|_ and e_ are orthonormal, one immediately verifies that 

L + R^R ~ e+ef e^e^ + e_ef e^e^ = e+e^ + e_e^ = 1 

almost surely. Let E e 5 be given, with 

s=(? ;_). (97) 

One finds that 

T(S) = E{LY.L^ + RY.R^\ = (e^I]e+ -f- e^I]e_)E[e,ef ] = ^-^-i-^E*, 
with the fixed point E.^ of T given by 

Ej; = 2 



q 
\-q 
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In particular, T(T,) = for traceless matrices S, i.e., for (t+ + (t_ = 0. Thus, Assumption [3] 
holds with K = 0. Define a weight function a; : — M by setting a;(Ae±) = A"* for all A € M, 
and uj{^) = for all ^ € that are not parallel to e+ or e_, with a value a; > 1 yet to be 

determined. For arbitrary ^ = (^i, ^2) G li^, it follows that 

5(0 E[w(L^O + ^{R^O] = 9(^(66+) + c^(ae-)) + (1 - g)(c^(6e+) + u;(6e_)) 

Now if C = Ae±, then = C| = AV4, and so 

9(0 = - l^iO- 
If instead ^ is not parallel to 6+ or e_ , then we estimate 

+ (l--g )C2i^i4 / 2max(g, 1 - q) 

Choosing cj 4max(g, 1 — g), it then follows that 



for all ^ E R'^, verifying Assumption [2] with p — 4 and Kp — 1/2. 

Remark 3. We emphasize that ([TU)) would not /jo/c? mi/i some Kp < 1 if we had made the 
canonical choice u!{£,) = |^|*. 

Consequently, the equation ([TJ possesses a unique stationary distribution /ioo of unit tempera- 
ture: fioo is supported on the coordinate cross and satisfies 

^iv)dn^{v) - (4^)-i/2 /" (g<^(u>, 0) + (1 - qMO, w))e-'""/Mw (98) 



for all test functions (p E (M^ ) . The fact that (IM)) defines an equilibrium is most easily checked 
on the level of its Fourier transform, 

*(»y) := Aoo(»/) = ge"''i + (1 - q)e-'''^. 
Indeed, ^E* satisfies the stationary equation (1^51) for every ^ S M^: 
E[<f{L^O^{R^^)] = gvl/(ae+)vl/(ae_) + (1 - q)^{i2e+)^{^2e^) 

= q{qe-^'i'-' + (1 - q)e-^'-'^f + (1 - q){qe-^i'^ + (1 - g)e-«^'/2)' = *(0. 

7. Appendix 

Here we list a variety of inequalities that are used throughout our calculations. 

• Young's inequality: for every p > 2, and for all x,y >0, 

^Pl^-^y < ^a;P/2 + lyPl^. (99) 

p p 

• Inequality on moments: for every s > 1, there is a constant Cs such that for all x, y > 0, 

{x + yY < + j/^ + C,{x'~^y + xy^-^). (100) 

• Remainder estimate for Fourier transforms: for p G (2, 3), and for all u e M, 

|e-™ - (1 - m)| < min{2|w|, < 2^~p\u\p/^, (101) 

|e-™ - (1 - m + < niin{M2, l\u\^} < 6^-p\u\p. (102) 
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• Von Bahr and Esseen inequality: Let Zi , . . . , Z„ independent (real valued) random 
variables such that E[Zi] — and E[|Zi|P] < +00 for some 1 < p <2, then 

n n 

E[|^Z,n < 2^E[|Z,|P]. (103) 

i=i 1=1 

See [22] for a proof. 
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